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Abstract 

Linear  plants  with  control  signals  entering  the  open  loop  dynamics  as  amplitudes 
of  periodic,  zero  average  functions  are  considered.  Stabilizability  properties  of  such 
plants  by  state  space  and  output  feedback  controllers  are  analyzed. 
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1.  INTRODUCTION  AND  PROBLEM  STATEMENT 


Consider  the  following  problem:  Given  the  system 

x  =  Ax  +  Bu f(t)  , 

x  €  IR\  u6  Rm ,  y  €  2RP  ,  (1.1) 

/(*)€«,  m  =  /(t  +  n  o, 

i  j*  f(t)dt  =  0  , 


find  a  state  space  controller  I\  £  IRmxn  such  that  with 

u  =  Kx  , 

(1.2) 

the  closed  loop  system 

x  =  (A  +  BKf(t))x 

(1.3) 

has  desired  dynamic  properties. 

A  characteristic  feature  of  this  problem  is  that  the  control,  u,  enters  the  open  loop 
dynamics  as  an  amplitude  of  a  periodic,  zero  average  function,  and  this  amplitude  can  be 
chosen  as  a  function  of  the  state  or,  more  generally,  output.  Systems  of  this  form  arise  in 
a  number  of  practical  applications.  One  of  these  is  the  Higher  Harmonic  Control  (HHC) 
in  helicopters.  Here  a  periodic  feathering  of  rotor  blades  around  a  fixed  pitch  angle  is 
introduced  in  order  to  suppress  the  fuselage  vibrations.  The  amplitude  of  the  vibrations  is 
typically  chosen  as  a  function  of  the  acceleration  of  pilot’s  seat  (system  output).  Recent 
experiments  [lj- [3]  have  shown  that  HHC  may  lead  to  an  order  of  magnitude  reduction  in 
fuselage  vibrations.  However,  no  formalized  methods  for  HHC  controllers  design  have  been 
described  in  the  literature. 

Another  example  is  the  periodic  operation  of  chemical  reactors  [4].  Here  again  the  prob¬ 
lem  !°  to  choose  the  amplitude  of  input  flow  vibrations  so  that  the  closed  loop  system  behaves 
as  desired. 
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From  the  theoretical  standpoint,  problem  (1.1),  (1.2)  is  closely  related  to  the  problem 
of  vibrational  control  [5]-[7]  which  can  be  formulated  as  follows:  Given  A  £  2Rnxn  find  a 
periodic,  zero  average  matrix  B(t)  £  IRn*n  such  that 

x  =  {A  +  B(t))x  (1.4) 

is  asymptotically  stable.  The  only  difference  is  that  (1.4)  does  not  contain  structural  con¬ 
straints  imposed  by  the  feedback  whereas  (1.3)  does.  Nevertheless,  due  to  the  obvious 
similarities  between  (1.3)  and  (1.4)  (including  the  methods  of  their  analysis-see  below),  the 
problem  (1.1),  (1.2)  will  be  referred  to  as  closed  loop  vibrational  control.  The  effect 
of  the  structural  contraints  on  the  problem  of  vibrational  stabilizability  is  the  topic  of  this 
note. 

The  results  presented  below  differ  also  from  those  of  [8]  in  that  the  latter  uses  the  vi¬ 
brations  introduced  in  the  parameters  of  dynamic  output  controllers  whereas  the  plants  are 
time  invariant. 

All  three  cases  however,  i.e.  [5]-[7],  [8]  and  the  present  work,  involve  linear  systems 
with  periodic  coefficients.  In  order  to  simplify  the  analysis  and  obtain  constructive  results, 
following  [5]-[8],  we  assume  that  the  periodic  function  f(t)  is  of  high  frequency  as  compared 
with  the  dynamics  of  x  =  Ax.  Formally,  this  means  that  function  /  has  the  asymptotic  form 
/(*)  where  0  <  t  «  1,  is  sufficiently  small.  Thus,  more  precisely  the  problem  addressed  in 
this  note  is  as  follows: 

Given 

x  =  Ax  +  Buf (-)  , 

y  =  Cx  ,  (1.5) 

x  £  lRn,  u  £  1R,  y  €  IR,  f  £  JR  , 

/(•)  periodic,  zero  average,  0  <  e  <<  1  small  parameter,  determine  under  what  conditions 
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there  exists  a  time  invariant  state  space  controller 

u  =  Kx  ,  (1.6) 

or  a  time  invariant  output  controller 

u  =  Kx  , 

x  =  Ax  +  Buf  +  L{y  -y)  ,  (1.7) 

y  =  Cx , 

such  that  the  closed  loop  dynamics  are  asymptotically  stable.  The  feedback  laws  (1.6), 
(1.7)  are  restricted  to  be  time  invariant  for  reasons  of  practical  implementation  and  possible 
uncertainty  in  the  knowledge  of  /(•).  Problem  (1.5),  (1.6)  is  considered  in  Section  2  and 
problem  (1.5),  (1.7)  is  discussed  in  Section  3.  In  addition,  we  characterize  the  pole  placement 
capabilities  ensured  by  closed  loop  vibrational  control  and  present  the  corresponding  results 
in  Section  4.  The  conclusions  are  formulated  in  Section  5. 

2.  STATE  SPACE  FEEDBACK 

Theorem  2.1:  There  exists  c0  >  0  such  that  for  all  0  <  e  <  to  system  (1.5)  is 
stabilizable  by  a  state  space  feedback  (1.6)  if  and  only  if  ( A,  B )  is  stabilizable  and  the  sum 
of  all  the  controllable  eigenvalues  of  A  is  negative. 

Proof:  Necessity  is  proved  by  the  following  considerations.  The  state  model  in  (1.5) 
has  the  Kalman  controllable  form 


Since  Anc  is  not  affected  by  feedback,  the  stabilizability  of  [A,B)  is  n  "'-css  ary. 

Consider  the  controllable  part  of  (2.1): 

xc  =  Acxc  +  B\uf  ,  xc  €  ZRm,  u  e  IR  .  (2.2) 
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Introducing  a  state  feedback  u  =  Kxc,  we  obtain 

Xe  =  (ac  +  BxKf  0j  xc  .  (2.3) 

Since  (2.3)  is  periodic,  there  exists  a  Lyapunov  transformation  which  reduces  (2.3)  to  an 
equation  with  constant  coefficients, 

i  =  Az  , 

preserving  the  stability  property.  The  following  equality  is  true  [5]: 

^  £  Tr  0  +  BxKf  0]  dt  =  Tr  A 
where  T  is  the  period  of  f(t/e).  Thus 

Tr  Ac  =  Tr  A  , 

where  Tr  Ac  is  equal  to  the  sum  of  all  the  controllable  eigenvalues.  This  completes  the  proof 
of  necessity. 


The  suffiency  is  proved  as  follows:  Consider  the  Kalman  controllable  form  (2.1)  of  the 
system  (1.5)  where  all  the  eigenvalues  of  Anc  have  negative  real  parts.  Without  loss  of 
generality,  assume  that  (2.2)  is  in  the  controllable  canonical  form  i.e.: 

xc  =  Acxc  +  B\uf  ,  xc6  IRm ,  u  €  JR 


where 


r  ° 

l 

0  ■ 

'  0  ' 

o 

II 

0 

. . . 

1 

,  Bi  = 

0 

L 

-fll  . 

.  1  . 

and  a,  are  the  coefficients  of  the  characteristic  polynomial  of  matrix  Ac.  Apply  state  feedback 


where  &,•  ~  1,  i=2,  . . .,  m.  In  the  fast  time  r  =  t/e,  the  closed  loop  system  is 


dxc 

dr 


{(■Ac  -)-  B\  K\f{r))  xc  . 


(2.5) 
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Let  $(i/e)  be  a  fundamental  matrix  for  (lle)BKif(t/e).  Reducing  (2.5)  into  the  standard 
form  [9]  and  then  applying  the  averaging  principle  we  have  the  averaged  system’s  equation 


(2.6) 


where 


0  1  •••O' 

0  _  0  _  •••  1 
~am  -  k2km<t>2  -am^i  —  k2km^i4>2  •••  —  ax  . 


and 

*(D=f  ^ 

Let  denote  the  open  loop  eigenvalues  of  (2.2)  and  choose  the  closed-loop  eigen¬ 

values  in  the  following  manner: 


1  m 

A,  =  -X>+;imA<.  (2.7) 

m  i=i 

We  determine  the  coefficients  al  of  the  closed-loop  characteristic  equations  corresponding  to 
A1; . . . ,  Am.  The  state  feedback  gains  (2.4)  can  be  found  by 


,  (a.  ~  a.) 

'  "  h 


i  =  2, . . . ,  m  . 


(2.8) 


The  control  gains  (2.8)  guarantee  the  asymptotic  stability  of  the  averaged  system  (2.6). 
As  it  follows  [8],  if  (2.6)  is  asymptotically  stable,  there  exists  eo  >  0  such  that  for  all 
0  <  e  <  Co  equation  (2.2)  is  also  asymptotically  stable.  This  proves  the  sufficiency.  Q.E.D. 

3.  OUTPUT  FEEDBACK 


Theorem  3.1:  There  exists  c0  such  that  for  all  0  <  c  <  e0  system  (1.5)  is  stabilizable 
by  an  output  feedback  (1.7)  if  (A,  B,C)  is  controllable  and  observable  and  Tr  A  <  0.  The 
separation  principle  holds,  i.e.  the  choice  of  K  and  L  can  be  carried  out  independently. 
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Proof:  Consider  system  (1.5)  with  {A,B,  C)  controllable  and  observable,  the  ob¬ 


server, 

x  =  Ax  +  Buf  +  L(y  —  Cx)  , 

and  the  feedback  law 


e 


In  fast  time  r  =  </c,  the  resulting  closed-loop  equations  are: 


dx 

dr 

'  tA  BKJ(t)  1 

X 

dx 

-dr  - 

tLC  e{A  -  LC)  +  BKxf(r)  \ 

X 

Let  $(f/e)  be  a  fundamental  matrix  for  (l/e)BKif(t/e).  Define 


(3.1) 


$ 


Reducing  (3.1)  into  the  standard  form  [9],  and  applying  the  averaging  principle,  we  obtain 
the  following  averaged  equations: 

’  *  1  _  A  +  $LC  -LC  -  A  +  LC  -  $LC 
x  ~  <&LC 


To  simplify  (3.2),  introduce  the  following  transformation 


X 

'  / 

0 

X 

e 

/ 

-I 

X 

which  yields 


x 

e 


— $1  A  —  LC  +  $LC 
0  A  -  LC 


x 

e 


(3-3) 


where  e  is  the  innovations  error.  Thus,  the  separation  principle  holds:  the  eigenvalues  of 
the  averaged  closed-loop  system  are  the  union  of  those  of  and  those  of  A  —  LC.  Using 
construction  similar  to  (2.6)-(2.8),  one  can  compute  the  state  feedback  gain  required  to 
stabilize  and  assign  all  the  eigenvalues  of  A  —  LC  through  the  choice  of  L  so  that  the 
stability  of  the  closed-loop  system  is  guaranteed.  Q.E.D. 
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4.  POLE  PLACEMENT  CAPABILITIES 


Consider  again  system  (1.5)  with  feedback  (1.6)  and  assume  that 

Tr  f 

A  =  -  =  |  —  ...  —  , 

e  [  e  c 

where  ki  ~  1,  t  =  1, . . . ,  n.  Thus,  the  closed  loop  system  is 

i  =  (A  +  BTf0)x ■ 

Using  the  averaging  theory  [9],  this  equation  can  be  reduced  to  the  averaged  equation, 

x  =  (A  +  B)x,  (4.2) 

and,  as  it  follows  from  [9]  and  [10],  (4.1)  is  asymptotically  stable  if  (4.2)  is  asymptotically 
stable.  The  stability  properties  of  (4.2)  can  be  checked  using  the  following: 

Theorem  4.1:  Assume  A  and  B  are  in  the  controller  canonical  form.  Then  there 


exists  c0  such  that  for  all  0  <  c  <  c0, 


0  •••  0  0 

0_  •••  0_  0 
—  bn  ■  •  ■  —  62  0 


00  lr 

T~  2  A,2  'W 


where  a,,  s  =  1,2,...,  are  the  Fourier  coefficients  of  /(t),  i.e. 

OO 

/(r)  =  ^o,sin(sr  +  v?s)  . 


Proof:  Follows  directly  from  Theorem  3  of  [5]. 


Q.E.D. 


It  follows  from  Theorem  4.1  that  the  characteristic  polynomials  of  A  and  {A  +  B)  have 
respectively,  the  form 


Po(s)  —  3n  -f-  a\oSn  1  +  d20Sn  2  +  •  ■  •  +  flnO  ) 
p(s)  =  Sn  4-  diSn  1  +  d2Sn  2  +  •  •  •  +  dn  , 


where 


ai  —  uio  ,  (4.6) 

a2  &20  1  (4.7) 


and  a},  3  <  j  <  n  can  be  arbitrarily  assigned.  Below  we  analyze  to  what  extent  the 
constraints  (4.6),  (4.7)  prevent  the  control  designer  from  assigning  the  closed  loop  eigenvalues 
in  a  desired  region  of  the  complex  plane.  More  specifically,  considering  the  closed  region 
D(cr,  <?)  of  Figure  4.1,  our  purpose  is  to  identify  the  conditions  under  which  we  can  find  n 
complex  numbers  occurring  as  pairs  of  complex  poles,  such  that  they  are  the 

roots  of  the  polynomial  (4.5)  which  satisfies  (4.6),  (4.7).  When  this  is  possible,  we  say  that 
pole  assignment  to  the  region  D(a.o)  using  vibrational  feedback  is  possible. 

As  preliminary  results,  we  have  the  following 

Lemma  4.1:  The  pole  assignment  in  the  region  D(a,^>)  using  closed  loop  vibrational 
control  is  possible,  only  if 

—  Ctio  <  TUT  .  (4.8) 


Proof:  If  Aj, . . . ,  An  €  D(a,  <t>)  are  the  .oots  of  the  polynomial  (4.5),  then 

-  a,  =  ^  A,  , 


i=i 


and 


Re(At)  <  a  . 


(4.9) 


(4.10) 


Equation  (4.8)  is  obtained  by  equating  the  real  parts  in  (4.9)  and  using  (4.10).  Q.E.D. 
Lemma  4.2:  Consider  the  optimization  problem:  Find  vector  r  6  IRm  such  that 

p(r)  =  H  r>rj  (4-11) 
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is  maximized  subject  to 


Z  r.  =  a  . 


1=1 


A  global  solution  r"  to  this  problem  is  given  by 


r  =  r2  =  .  .  .  =  rm  =  —  , 

m 


and 


(4.12) 


(4.13) 


(4.14) 


Proof:  S  inre  the  constraint  (4.12)  is  ab/ays  regular,  we  apply  the  Lagrange  multiplier 
rule.  The  Lagrangrian  is 

*n  /  m  \ 

L(r,  \)  =  £  r.r,  +  A  (  ]Tr,  -  aj  (4.15) 

•<j 

The  first  order  necessary  conditions 


dL{r,  A) 
dr 

dL(r, \) 
d\ 


(4.16) 

(4.17) 


yield  a  linear  system  of  equations  whose  solution  is  (4.13)  together  with: 


1  -m 

A  = - a  . 

m 


(4.18) 


The  second  order  conditions  ensure  that  (4.13)  is  a  strict  maximum.  Moreover,  since  un¬ 
der  the  constraint  (4.12)  the  cost  function  (4.11)  is  quadratic,  this  maximum  is  the  globa' 
maximum,  and  the  proof  is  complete.  Q.E.D. 


We  can  now  state  the  main  result: 


Theorem  4.2:  The  real  pole  assignment  in  the  region  D(cr,(f>)  i°  possible  using 
vibrational  feedback  if  and  only  if 


—  ciio  <  i  (4-19) 

«20  <  “T“aio  •  (4-2°) 
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Proof:  The  necessity  of  (4.19)  follows  from  Lemma  1.1.  The  necessity  of  (4.20) 
follows  from  the  fact  that  when  rt,....rn  are  the  roots  of  (4.5),  the  coefficients  ax  and  a2 
are 

n 

«i  =  ~Er'  ’  (4.2!) 

1  =  1 
n 

«2  =  E  W  -  (4.22) 

*,J  =  1 
«>J 

The  maximum  value  that  (4.22)  can  achieve  subject  to  the  constraint  (4.21),  (4.6)  given 

by  Lemma  4.2  and  is  exactly  the  right  hand  side  of  (4.20).  Therefore  ;f  (4.20)  is  violated, 

pole  assignment  to  the  region  D(a,d >)  with  real  poles  is  not  possible. 


To  prove  the  sufficiency  of  (4.19),  (4.20),  assume  they  both  hold.  Choose 


aio 

n  =  r2  =  . . .  =  rn  =  —  . 

n 


(4.23) 


It  L  .mmediately  checked  that  these  real  numbers  solve  the  problem  of  pole  assignment  to 
the  region  D(cr,<t>),  which  completes  the  proof.  Q.E.D. 


Theorem  4.2  gives  a  simple  solution  of  the  problem  of  pole  assignment  to  the  region 
D(cr,  d>)  with  real  poles.  When  non  real  poles  are  allowed,  more  complicated  results  similar 
to  Theorem  4.2  can  be  derived,  and  will  be  included  in  a  more  entensive  version  of  this 
paper. 

5.  A  CONCLUDING  REMARK 


As  it  follows  from  the  above  and  [5],  the  conditions  of  closed  loop  vibrational  stabiliz- 
ability  are  remarkably  similar  to  the  conditions  of  open  loop  vibrational  stabilizability:  the 
only  difference  is  tho.t  the  former  requires  the  stabilizability  property  of  ( A,  B ).  Roughly 
speaking,  the  reason  for  this  is  that  a  controllable  pair  can  always  be  transformed  into  the 
controllable  canonical  form  and  this  transformation  removes  the  constraint  on  the  structure 
of  the  input  matrix  B. 
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